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ABSTRACT: 

This paper presents iterative determinant method for solving eigenvalue problems. A matlab program that 
operates iterative to evaluate the determinant of tlie problem was written. In the program, a trial 
eigenvalue is used in the program to compute the determinant. A small value (iterator) is added to trial 
eigenvalue to obtain a new trial eigenvalue, which is used to compute new determinant. The trial 
eigenvalue in the sequence that made the determinants to move from negative values to positive value or 
move from positive values to negative value becomes required eigenvalue. Some engineering problems 
were used to validate the eigensolver. Some of the data from the eigensolver and the corresponding exact 
data are 2.468 and 2.4678; 9.876 and 9.875; 60.001 and 60.00; 12.37 and 12.3695. Close look at 
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I. INTRODUCTION 

According to Ibearugbulem et al (2013), the stiffness matrix [k], the geometric matrix [kg] and the 
mass (inertia) matrix [ki] are formulated using the assumed shape function. This shape function is usually 
assumed to approximate the deformed shape of the continuum. If the shape function assumed is the exact one, it 
means the solution will converge to the exact solution. The inertia matrix and the geometric matrix formulated 
using the assumed shape functions are called consistent mass matrix and consistent geometric matrix 
respectively (Paz, 1980 and Geradin, 1980). In the work of Ibearugbulem et al (2013), the dynamic equation in 
structural dynamic as: 

[k]-Mki] = o (l) 

For static stability problem, the eigenvalue equation is given as: 
[k]-Nc[kg] = 0 ~ (2) 

Where K is the material stiffness matrix, k; is the matrix of inertia, k g is the geometric matrix X is natural 
frequency parameter and N c is the critical buckling load parameter. Of note here is that both k g and k t are 
consistent matrices (that is they are non-diagonal matrices). The difficulty posed by this type of eigenvalue 
problem led many researchers to transform the consistent matrices to diagonal matrices as: 
M-U[I] = 0 ' (3) 
[k] - Nc A[I] = 0 (4) 

Here, A[I] is the transformed diagonal inertia or geometric matrix. In dynamics, the diagonal inertia 
matrix is often called lumped mass matrix. Many methods are adopted so far by researchers for solutions of 
eigenvalue problems of equation (3) and (4). According to Ibearugbulem et al, (2013), the methods include 
Jacobi method, Polynomial method, Iterative method and Householder's method were used by (Greenstadt, 
1960; Ortega, 1967; and James, Smith and Wolford, 1977. Others are Power method, Inverse iterative 
(Wilkinson, 1965), Lanczos method (lanczos, 1950), Arnoldi method (Arnoldi, 1951; Demmel, 1997; Bai et al, 
2000;Chatelin, 1993; and Trefethenand Bau, 1997), Davidson method, Jacobi -Davidson method (Hochstenback 
and Notay, 2004; and Sleijpen and Vander Vorst, 1996), Minimum residual method, generalized minimum 
residual method were used by Barrett et al, (1994), Multilevel preconditioned iterative eigensolvers (Arbenz and 
Geus, 2005), Block inverse-free preconditioned Krylov subspace method (Quillen and Ye, 2010), Inner-outer 
iterative method (Freitag, 2007), and adaptive inverse iteration method (Chen, Xu and Zou, 2010), Matrix 
iterative-inversion (Ibearugbulem et al, 2013). Sadly, of all these methods, only polynomial and iterative- 
inversion methods can handle the problems of equations (1) and (2). Other methods can only be used for 
equation (3) and (4). However, polynomial method also becomes very difficult to use when the size of the 
matrix exceeds 3x3. In the same way, iterative-inversion method is also difficult when the matrix size is large 
and the speed of computation is very slow. 

The essence of this paper is to present a method that can be used in solving eigenvalue problems of equations 
(1), (2), (3), and (4) for any size of matrix with high speed and good accuracy 
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IT. ITERATIVE DETERMINANT 

In this method, a trial eigenvalue X 0 or Nc 0 of default value of zero shall be substituted into the 
eigenvalue equation to obtain eigenvalue matrix as: 
[k]- X[ki] = [kk] 0 (5) 

The determinant of the eigenvalue matrix, [kk] 0 shall then be computed and the value kept. The value 
of this initial determinant, dtO may be positive or negative. Now the value of the default eigenvalue shall be 
increased by adding iterator (say 1, 0.1, 0.01, 0.001 etc. as desired) to it to obtain new eigenvalue, X\. This shall 
be substituted into the eigenvalue equation to obtain [kk]j. The determinant, dtj of [kk]j shall be computed. If 
dto is negative and dtl is also negative, then the actual eigenvalue has not been obtained and the iterator has to 
be added to Xj to obtain X 2 - This process has to be continued until we reach a stage where X n . { > 0 and A. n < 0 or 
^ n _i < 0 and X n > 0. Note here that the bigger the iterator, the faster the computation and less the accuracy, the 
smaller the iterator, the slower the computation and more the accuracy. To ease the use of this method, a general 
matlab program was written. The user is at will to modify the iterator, f and sub iterators ff (i) in the program, 
where is 1, 2, 3 ... n and n is the size of the square matrix involved. For instance, if the order of the expected 
eigenvalue is 100, iterator and sub iterators of 1 or 0.1 can be used. In this case the expected accuracy shall be of 
the order of iterator, 1 or 0.1. In the same way, if the order of expected eigenvalue is 1, iterator of 0.001 can be 
used. The choice of iterator is guided by the speed and accuracy of the computation. In all, for high accuracy, 
iterators of 0.001 and below are advisable but the speed may be very slow. Thus, the level of accuracy and speed 
shall guide you in using your choice iterator in this program. The user is also at liberty to choose the range of 
eigenvalues to compute for by adjusting "for r = 1:3" to say "for r = 1:5" if the size of the matrix is 5x5. The 
stiffness matrix "prs" and the geometric matrix or inertia matrix "pri" in the program are to be entered by the 
user. The iterating count ceiling, "k" in the program can also be raised from 6260 to any higher value to suit the 
user. Number of sub iterators "ff(i) as we have them in the program is 7. The user is at liberty also to introduce 
more sub iterators say ff(8), ff(9) etc. as matches the size of the matrix. Some engineering problems were used 
to ascertain the adequacy of the method and the program. 



III. NUMERICAL EIGENVALUE PROBLEMS 

The program will be used to test the following problems. 
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IV. RESULT AND DISCUSSION 

Table 1 shows result data from eigenvalue problems herein. Eigenvalues obtained from the program were 
compared with the exact eigenvalues. The exact eigenvalues were obtained by trial and error means using 
Microsoft excel worksheet. Any eigenvalue that made the determinant of the eigenvalue matrix [kk] to become 
zero is the exact eigenvalue. The data from the program compared very well with the values from Microsoft 
excel worksheet with high degree of accuracy. As said earlier, it was observed that with small iterators say 0.01, 
we obtained more accurate eigenvalues with slow computing speed. We also confirmed that with big iterators 
say 1.0 we obtained less accurate eigenvalues with fast computing speed. 



Table 1 : Result data from the Eigenvalue Problems 
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V. APPENDIX A (MATLAB PROGRAM) 

nn = input('enter size of matrix');nn = nn*l; 

p=0;m=l ;k=6260;j=l ;f=0. 1 ;ff(l)=0. l;ff(2)=0.01 ;ff(3)=0.001;ff(4)=0.001 ;ff(5)=0.001 ;ff(6)=0.001 ;ff(7)=0.001; 
f or r = 1 : 3 
while p < k 

prs = [43.2 0 0 ;0 6.400093 0 ;0 0 6.4]; 

pri = [0.54824 0 0 ;0 0.01268 0 ;0 0 0.01268]; 

a = prs - p*pri; 

for x = l:nn 

for y = l:nn 

c(x,y)= a(x,y) ; 

end 

end 

d = det(c); 

if (m< 1.1); tl =d;end 

m = m + 1; 

if (j >nn); break; end 

if ((d >= 0)&& (tl <= 0));py(j) = p;j = j + l;tl = d;end 
if ((d <= 0)&& (tl >= 0));py(j) = p;j = j + l;tl = d; end 
p = p + f; 
end 

p = p-f; f=ff(r);m=l; 
end 
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